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Abstract 

This talk describes the reasons why t] and rj' decays are an interesting topic 
of study for both theory and experiment. The main part discusses the re- 
sults of the recent calculation of rj ^ at two-loop order in ChPT. Some 
puzzling aspects of the results compared to earlier dispersive calculations are 
highlighted. I also like to remind the reader of the use of rj and rj' decays for 
studying the anomaly. 
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Abstract 

This talk describes the reasons why rj and rj' decays are an interesting topic 
of study for both theory and experiment. The main part discusses the re- 
sults of the recent calculation of — Svr at two-loop order in ChPT. Some 
puzzling aspects of the results compared to earlier dispersive calculations are 
highlighted. I also like to remind the reader of the use of t] and t]' decays for 
studying the anomaly. 

1 Introduction 

This conference has a lot of talks related to rj and rj', both on decays, produc- 
tion and in a hadronic medium. The production is treated in a plenary talk 
by Krusche and decays experimentally in the talk by Wolke. There were also 
a lot of talks for both production and decay in the parallel sessions. In this 
talk I will concentrate on decays and in particular mainly on 77 — > Sir. There 
are lots of references treating r] and t]' physics. Many of them can be found in 
the proceedings of two recent conferences devoted to them [1,2]. There have 
also been more recent workshops in Jiilich (ETA06) and Peniscola (ETA07). 

This talk first discusses why rj and f]' are interesting, then reminds the 
reader of some of the aspects of Chiral Perturbation Theory (ChPT) after 
which the main part, devoted to rj ^ Sn comes. I close by pointing out 
some properties of rj' rj-K-K^ Tmn decays and the anomaly. Earlier reviews 
covering similar material are Refs. [3,4]. 

2 Why are 77 and rj' Interesting? 

The T] and t]' are particles that decay strongly but all their decays are sup- 
pressed. That means that they are good laboratories to study non-dominant 
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strong interaction effects. Weak decays can happen but do occur at branch- 
ing ratios of order 10^^^ or lower. So, if charge conjugation violation would 
be discovered it would be very important. On the other hand, most standard 
extensions of the standard model do not predict such effects at an observable 
level in r] or r]' decays. 

But let us first see why pseudo-scalars are special. The QCD Lagrangian 

is 

Cqcd = [-^L^ Ql + iqRip Qr - rriq {q^qi + Ql^r)] (1) 

q=u,d,s 

So if niq = then the left and right handed quarks are decoupled and they 
can be interchanged freely among themselves leading to a global symmetry 
G = U{3)lxU{3)r. This symmetry is clearly broken in the hadron spectrum, 
the proton and the Sn, as well as the p and the ai have very different massed]. 
The chiral symmetry group G must thus be spontaneously broken, only the 
vector part of the group is clearly visible in the spectrum. 

As a consequence there must be a set of light particles, the pseudo- 
Goldstone boson, whose interactions vanish at zero momentum as follows 
from Goldstone's theorem. There are eight fairly light particles around with 
the right quantum number, 7r°, vr^, K^, K^, and rj. But the next candi- 
date with the correct quantum numbers, the f]', is heavy. We write the group 
G in terms of simple groups, 

G = U{3)l X U{3)n = SU{3)l x SU{3)n x U{l)v x U{1)a , (2) 

and notice that the breaking pattern of G = SU{3)l x SU{3)r — > H = 
SU{3)v gives eight light particles as observed. The reason is that the U{1)a 
part of G is a good symmetry of the classical action but not of the full 
quantum theory. The divergence of its current has a part coming from the 
anomaly which couples to gluons via 

d^A'^ = 2^fUj with u = ^e^''''HiG^,G^f,. (3) 

LO is strongly interacting thus the divergence of the singlet axial-current can- 
not be treated as zero. So the 1]' can be heavy as is seen experimentally. 
Quantum effects break thus the U{1)a, however the r.h.s. of (E]) is a total di- 
vergence, so how can it have an effect? The answer was found by 't Hooft [7]. 
Gauge field configurations with non-zero winding number, instantons, can 
produce an effect. This in turn led to the so-called strong CP problem but 
solved the r]' mass problem. A conclusion is thus that the t]' has potentially 

"'^ There is a discussion at present whether chiral symmetry is restored for higher hadron 
masses. This is not relevant for this talk. Recent references can be traced from [5,6]. 
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large and very interesting non-perturbative effects and interactions with glu- 
onic degrees of freedom that differ from other hadrons. Since ifi ^ rus this 
also affects 77 physics via mixing. 



3 Chiral Perturbation Theory 

The chiral symmetry of QCD and its spontaneous breaking has many con- 
sequences. The best method to exploit these is Chiral Perturbation Theory 
(ChPT) which is best defined via 

ChPT = "Exploring the consequences of the chiral symmetry of QCD and 

its s])()iitaii(M)us l)reakiiig using ofToctivo field th(M)ry toc4nii(iuos." 



A derivation which clearly brings out all the assumptions involved is [8]. 
Lectures and review articles can be found in my LatticcOT talk [9] or on the 
webpage [10]. The original modern references are [11,12]. 

ChPT uses as power-counting essentially dimensional counting in terms of 
a generic momentum p. Momenta and meson masses are counted as order p. 
Because of the Gell-Mann-Oakes-Renner relation, m\j oc m^, quark masses 
and external scalar and pseudo-scalar fields are counted as order and 
the covariant derivative requires external vector and axial-vector field to be 
counted as order p. With this counting there is no term of order p° in the 
chiral Lagrangian. The lowest order Lagrangian is given by 

A = ^{{D.U'^D^U) + (x^C/ + xC/^)} , (4) 
with U parameterizing the Goldstone Boson manifold G/H with 



U{(t>) = exp(iv^$/Fo) , 



— F H F 7r+ A + 



V2 V6 



K- K 



2vs 



\ V6 / 

(5) 

The external fields are in the covariant derivative, -D^t/ = d^U — ir^U + iUl^ 
for the left and right external currents: r(/)^ = + (— )a/i and the external 
scalar and pseudo-scalar external densities are in x = 2i?o(s + ip). Quark 
masses come via the scalar density s — A4 + ■ ■ ■ and traces are over (quark) 
flavours {A) — Trp {A). The number of parameters increases fast at higher 
orders, there are 10+2 at order [13] and 90-1-4 at order p^ [14] for three- 
flavour mesonic ChPT. 
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Vtt In 
/TT \7r 

Figure 1: Resonance saturation of the order low-energy-constants C[ via 
resonate exchange for vrTr-scattering. 




The main uses of ChPT are that it contains all the SU{3)v relations 
automatically and in addition relates processes with different numbers of 
pseudo-scalars and it includes the nonanalytic dependences on masses and 
kinematical quantities, often referred to as chiral logarithms. As an example, 
the pion mass in two-flavour ChPT is given by [12] 



rrii = 2Bm + 



327r2 ° /i2 



+ •••, (6) 



with = 2Bm and B ^ Bq, F ^ Fo because of two versus three-flavour 
ChPT. In ([6]) we see the logarithm and the occurrence of the higher order 
parameter iKfi). Eq. also shows some of the choices that need to be 
made when performing higher order ChPT calculations: Which subtraction 
scale fi and which quantities should be used to express the results. Lowest 
order masses or physical meson masses and dito for the decay constants 
and other kinematical quantities as s, t, u in vrvr-scattering. There is clearly 
no unique choice and the choice can influence the apparent convergence of 
the ChPT series quite strongly. Another problem is that typically, not all 
the higher order parameters that show up in the calculations are known 
experimentally. Thus one needs to make estimates of these, usually via a 
version of resonance saturation originally introduced in ChPT in [15]. This 
is schematically depicted in Fig. [H More recent references on resonance 
saturation and possible pitfalls are [16,17]. Discussions on this problem can 
also be found in the papers on order p^ ChPT and the review [18]. 



4 77 ^ Stt 

In the limit of conserved isospin, i.e. we turn off electromagnetism and set 
m„ = m^, the rj is stable. Direct electromagnetic effects have been known 
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to be small since long ago [19,20]. It should thus proceed mainly through 
the quark-mass difference — rrid- The lowest order was done in [21,22], 
order in [23] and recently the full order has been evaluated [24] . In this 
section I will mainly present the new results of [24]. 

The momenta for the decay rj — >• 7r"'"7r~7r° are labeled as p.^, p+, p_ and 
Po respectively and we introduce the kinematical Mandelstam variables 

s = (P+ + P-Y ,t= {?++ Pof ,u = {p^+ PqY . 

These are linearly dependent, s + t + u = m^o + + + = Ssq .The 
amplitude is for the charged and neutral decay 

(vr^Tr+Tr^outl?]) = i {27r)'^ 6^ {pr, — PTr+ — P-k- — P-ko) A{s,t,u) , 
(7r°7r°7r°out|r7) = i {27i f 6^ {p^ - pi - p2 - Ps) A{si, S2, S3) , 

^(■51, 52,53) = A{si,S2,S3) + A{s2,S3,Si) + A{s3,Si,S2) . (7) 

The relation in the last line of (JTj) is only valid to first order in — m^. The 
factor of rriu — ma can be pulled out in various ways. Two common ones are 

A{s,t,u) = i-M{s,t,u) or A{s,t,u) = —^{ml - ml) ^''^ 



(8) 

with quark-mass ratios i? = {ms—m)/{md—mu) and = R{ms+'md)/{2rh). 
The lowest order result corresponds to 

M(s,t,«Uo= ((4/3)m2-,)/F2. (9) 

The tree level determination of R in terms of meson masses gives with a 
decay rate of 66 eV which should be compared with the experimental results 
of 295±17 eV [25]. In principle, since the decay rate is proportional to l/R"^ 
or 1/Q^, this should allow for a precise determination of R and Q. However, 
the change required seems somewhat large. The order calculation [23] 
increased the predicted decay rate to 150 eV albeit with a large error. About 
half of the enhancement in the amplitude came from vrvr rescattering and the 
other half from other effects like the chiral logarithms [23]. The rescattering 
effects have been studied at higher orders using dispersive methods in [26] 
and [27]. Both calculations found a similar enhancement in the decay rate 
bringing it to about 220 eV but differ in the way the Dalitz plot distributions 
look. This can be seen in Fig. [2] where I show the real part of the amplitude 
as a function of s along the line s = u. The calculations use a very different 
formalism but make similar approximations, they mainly differ in the way the 
subtraction constants are determined. That discrepancy and the facts that in 



5 



Johan Bijnens 



1] and rj' physics 




Figure 2: (a) Decay amplitude obtained by use of extended Khuri-Treiman 
equations [26] along the line s = u. (b) Alternative dispersive analysis for 
the decay amplitude [27]. Figs, from [4], adapted from [26,27]. 

Km the dispersive estimate [28] was about half the full ChPT calculation [29] 
and at order p'* the dispersive effect was about half of the correction for 
?7 — i> Svr makes it clear that also for this process a full order calculation is 
desirable. This has been done recently in [24]. 

Ref. [24] generalizes the methods of [30] to deal with tc^-i] mixing to 
processes with mixing on more than one external leg. The input parameters 
are from the main order fit, called fit 10, of [30] and the needed order p^ 
constants are determined by resonance exchange as discussed earlier. Details 
can be found in [24]. In Fig. [S]! show the numerical result for the amplitude 
along two lines in the Dalitz plot, t = u and s = u. The latter can be 
compared directly with the dispersive result of Fig. O The correction found 
in [24] at order p^ is 20-30% in amplitude, larger in magnitude than the 
dispersive estimates [26,27] but with a shape similar to [27]. 

The Dalitz plot in 77 — 37r is parameterized in terms of x and y defined 
in terms of the kinetic energies of the pions Tj and Qn = — 2m^+ — m^o 
for the charged decay and z defined in terms of the pion energies Ei. The 
amplitudes are expanded in x,y, z. 

X = ./3 ^+ ~ ^ , i/ = ^-i, z=^-y (^^^^^\\ 

Qrj Q 

|M(s,t,M)p = Al{l + ay + by'^ + dx^ + fy"^ + ■■■) , 

|M(s,t,M)p = Zo(l + 2«2 + ---) • (10) 
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L: fit 10 and C: 
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Figure 3: (a) The amphtude M{s, t, u) along the hne t = u. The vertical 
lines indicate the physical region. Shown are the real and imaginary parts 
with all parts summed up to the given order, (b) Similar plot but along the 
line s = u. Figs, from [24]. 



Table 1: Measurements of the Dalitz plot distributions in 77 — tt+tt 7r°. The 
KLOE result [31] for / is / = 0.14 ± 0.01 ± 0.02. 



Exp. 



KLOE [31] 
CB [32] 



[34] 



-1.090±0.005t°:[!°^ 
-1.22 ±0.07 
-1.08 ±0.014 
-1.17±0.02 



0.124±0.006 ± 0.010 0.057±0.006t[!:°°^ 

0.22 ± 0.11 0.06±0.04 (input) 
0.034 ±0.027 0.046 ±0.031 

0.21 ±0.03 0.06 ±0.04 



Recent experimental results for these parameters are shown in Tabs. [T] and 
[21 There are discrepancies among the experiments but the two latest preci- 
sion experimental measurements of a agree. The predictions from ChPT 
to order with the input parameters fixed as described earlier are give in 
Tabs. [3] and |H The predictions from the dispersive analysis as well as [38] 
have not been included. The different lines corresponds to variations on the 
input and the order of ChPT. The lines labeled NNLO are the central re- 
sults. The agreement with experiment is not too good and clearly needs 
further study. Especially puzzling is the a is consistently positive while the 
dispersive calculations as well as [38] give a negative value. The inequality 
a < {d + b — a'^/4:) /4 derived in [24] shows that a has rather large cancel- 
lations inherent in its prediction and that the overestimate of 6 is a likely 



7 



Johan Bijnens 



1] and 7]' physics 



Table 2: Measurements of the Dalitz plot distribution in 77 — > 7r°7r°7r°. 



Exp. 


a 


KLOE [35] 
Crystal Ball [36] 
WAS A/CELSIUS [37] 


-0.027 ± 0.0041° °°^ 
-0.031 ± 0.004 
-0.026 ±0.010 ±0.010 



Table 3: Theoretical estimate of the Dalitz plot distributions in 77 ^ tt+ tt 7r°. 





42 


a 


b 


d 


f 


LO 


120 


-1.039 


0.270 


0.000 


0.000 


NLO 


314 


-1.371 


0.452 


0.053 


0.027 


NLO m = 0) 


235 


-1.263 


0.407 


0.050 


0.015 


NNLO 


538 


-1.271 


0.394 


0.055 


0.025 


NNLO (/X = 0.6 GeV) 


543 


-1.300 


0.415 


0.055 


0.024 


NNLO In = 0.9 GeV) 


548 


-1.241 


0.374 


0.054 


0.025 


NNLO (C[ = 0) 


465 


-1.297 


0.404 


0.058 


0.032 


NNLO m = C[ = 0) 


251 


-1.241 


0.424 


0.050 


0.007 



cause of the wrong sign obtained for a. In addition, the fairly large cor- 
rection obtained gives in the end somewhat larger values of Q compared to 
those derived from the masses [24] . 

5 Other Remarks 

I would simply like to repeat here some remarks made earlier, see e.g. [4]. The 
hadronic decays of the 1]' are interesting, they are predicted to be small at 
lowest order, rj' Sir agrees reasonably well with expectations but rj' — > t^tttt 
has very large higher order corrections since the lowest order is suppressed 
by a factor of m^. I would also like to emphasize once more that the decay 

Table 4: Theoretical estimates of the Dalitz plot distribution in 77 — > Tr^TT^Tr". 





^2 
Aq a 


LO 
NLO 

NLO m = 0) 

NNLO 
NNLO (C[ = 0) 
NNLO m = C[ = 0) 


1090 0.000 
2810 0.013 
2100 0.016 
4790 0.013 
4140 0.011 
2220 0.016 
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of 1] and rj' allow many tests of the triangle, quadrangle, . . . anomaly. 
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